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ABSTRACT 

After recent AIA observations by Savage, McKenzie and Reeves we revisit the scenario proposed by 
us in previous papers. We have shown that sunward, generally dark plasma features originated above 
posteruption flare arcades are consistent with a scenario where plasma voids (which we identify as supra- 
arcade reconnection outflows, SAROs) generate the bouncing and interfering of shocks and expansion 
waves upstream of an initial localized deposition of energy which is collimated in the magnetic field di- 
rection. In this paper we analyze the multiple production and interaction of SAROs and their individual 
structure that make them relatively stable features while moving. We compare our results with observa- 
tions and with the scenarios proposed by other authors. 

Subject headings: 



1. Introduction 

Sunward dark moving trails with origin [40 - 
60] Mm above posteruption flare arcades and decel- 
erating speed in the range ~ [50 - 500]km s" 1 were 
first detected with the Yohkoh Soft X-ray Telescope 
(SXT). Since then, they have been extensiv ely reported 



using other instruments suc h as TRACE (llnnes et al 



l2003allbb . SOHO/SUMER (llnnes et aD l2003bh an d 
SDO/AIA dSavage et alJ l2012h . dWarren et al.l 1201 lb . 
The lack of X-ray and extreme-ultraviolet (EUV) 
signatures in images and spectra has lead to consen- 
sus that these downward moving structures are likely 
voided flows generated by reconnection processes in 
a current sheet above the flare arcade. Besides the 
dark moving structures, bright supra-arcade down- 
flowing features have also been reported during flares 
dMcKenziJ200Cb . 
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McKenzie and Savage! 12009b J Savage and McKenzie 
d201ll) suggested that supra-arcade downflows (S ADs) 
are the cross-sections of thin and empty flux tubes re- 
tracting from a reconnection site high in the corona. 
According to these authors, the high enough inner 
magnetic pressure could be the reason the voids are 
able to resist being fille d in immediately b y the sur- 
rounding denser plasma. iLinton et aD ( 2009|) proposed 
a scenario where the dynamic of retracting magnetic 
fields is triggered by a localized reconnection event 
that produces up and down flowing reconnected flux 
tubes, which are slowed down by underlying magnetic 
arcade loops. The observed SAD speeds are lower 
than expected for reconnection outflows in regions 
of typically assume d Alfven speeds of 1000 km s" 1 . 
Linton et al .] d2009h suggested that drag forces could 



work against the reconnection outflow. 

Verwichte etail d2005k analyzed TRACE downflow 
oscillations transverse to the magnetic field. They 
found that the initial speeds and the displacement am- 
plitudes, of kink-like type in the observational dark 
lanes of variable sizes (between [~ 2 - 9] Mm), de- 
crease as they propagate downwards, while the period 
remains fairly constant with height. 

Recently, after AIA observat i ons w ith high resolu- 
tion and cadence, ISavage et al.1 (|2012|) re-interpreted 
SADs as wakes created by the retraction of thin loops 
(SADLs) instead of the previous interpretation as flux 
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tube cross-sections. SADLs are features of sizes 
[~ 0.9 - 1.3]Mm observed during the early phase 
of the eruptive event and SADs as features of sizes 
~ 9Mm that become apparent afterwards, in contrast 
to bright, high-temperature plasma associated with 
current sheets. They proposed that deceleration is ex- 
pected due to the buildup of downstream magnetic 
pressure and/or drag mechanisms. 

In [Costaet al.1 (|2Q09|) , ISchulz et al.1 (l2010h and 
iMaglione et al.1 d201lh (hereafter Paper 1-3, respec- 
tively) we showed, by means of 1 and 2D MHD simu- 
lations, that the dark tracks are consistent with plasma 
voids generated by the bouncing and interfering of 
shocks and expansion waves upstream of the initial lo- 
calized deposition of energy. The composition of both, 
a resulting sunward directed hydrodynamic shock pat- 
tern and a perpendicular magnetic shock, produce an 
overall transversely shaking void that propagates to- 
wards the surface o f the sun, resembling th e kink-like 



mode described in Verwichte et al. (2005). Contrary 



to the ID results -where the sunwards dynamic is in- 
dependent of the magnetic field intensity owing to its 
exclusive waveguide role- in the 2D simulation the 
sunwards speed is higher for higher values of the mag- 
netic field. This can be interpreted as the capability of 
the low coronal plasma to collimate the deposition of 
energy in the magnetic field direction. 

As the obtained J3 values were larger than those out- 
side the voided cavity (see Fig. 2 of Paper 2), with the 
contour in total pressure equilibrium, we concluded 
that the internal magnetic pressure cannot be respon- 
sible of preventing the collapse of the vacuum zone. 
Instead, we found that the higher values of J3 inside the 
dark lanes are due to hotter plasma -with low inner val- 
ues of density (~ 2 10 7 cm" 3 )- conc entrated in regions 
of T > 20 MK IWarren et al.lll999h . Actual observa- 
tions do not show significant emission values of the 
voided plasma in any of the instrumental wavelengths 
(e.g., high-temperature AIA bandpasses). Hence, the 
observations could be consistent with these plasma pa- 
rameters. 

According to ISavage etaD d2012|) shrinking loops 
or SADLs are always leading SAD features. Both 
SADs and SADLs are interpreted as the outflows cre- 
ated during the re-organization of the magnetic field 
during the reconnection process. However, descend- 
ing SADs are only seen trailing loops once the cur- 
rent sheet has become noticeably filled with plasma. 
In this paper, in light of the new AIA observations 
we revisit the scenario proposed in Papers 1-3. We 



simply reinterpret our previous ID and 2D resulting 
voided features as supra-arcade reconnection outflows 
(SAROs) to diff erentiate them from the SAD descrip- 
tion as stated by lSavage et al.1 (|2012|) . In our 2D model 
we are not simulating the arcade loops nor the possible 
interaction of the SAROs with the loops, and will come 
back to the SAD and SADL features later. Hence, re- 
connection is modelled as a pressure pulse represent- 
ing the energy deposition that triggers the phenomenon 
that creates the dark track which we call SAROs. 

With the aim of testing the scopes of our model, we 
first focus on reproducing the dynamic and oscillat- 
ing interaction between several of these SARO voided 
structures. In Paper 1-3 we found that to reproduce the 
oscillating pattern we required closed boundary con- 
ditions. However, in Paper 2 we found that the os- 
cillating pattern was robust under a broad variation 
of the boundary conditions -contemplating different 
density values at the extremes and allowing total or 
partial boundary rebounds. Thus, the question arises 
about the role played by the boundary conditions. The 
observations show that the outflows sustain their in- 
dividual configuration all the way down into the ar- 
cade structure. Thus, the simulation of multiple out- 
flows triggered by different pressure pulses must main- 
tain these individualities while they evolve decelerat- 
ing sunwardly. However, can we consider the dynamic 
of a multiple outflow structure as the superposition of 
individual ones (as described in Paper 1-3)? If this is 
the case, what p r oduce s the wavy features described in 
Verwichte et al.1 d2005h? 



2. Numerical code and initial conditions 

We carried out the 2D-MHD simulations with the 
Mezcal code, which has been e xtensively tested in 
severa l sce narios (see for i nstanc e De Colle and Ragal 
d2006h or be Colle et aD d2008k ). All calculations 
were performed with a numerical grid of (x,y) = 
(400, 800) grid-points and a physical size of (20, 40)Mm, 
corresponding to a resolution of (50,50)km. The 
coordinate y represents the sunwards direction and 
the x coordinate the transverse to the magnetic field 
one. We assumed a constant radial initial magnetic 
field structure, and typical background temperature 
of T = 1.0 10 6 K. Simultaneous (t = Osec) spherical 
pulses (AP/P)i = (110,90,110) i = 1,3 (AP is the 
triggering pressure pulse and P is the background gas 
pressure of the corona) of a radius 0.6Mm were local- 
ized in positions (67,720), (167,680) and (333,720) 
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corresponding to (3.35, 36.0)Mm, (8.35, 34.0)Mm and 
(16.65, 36.0)Mm, respectively. 

Several simulations were carried out, varying the 
magnetic field, B and the density p. Also, we added 
a localized deposition of energy -modeled as a new 
triggering pressure pulse, (AP/P) 4 - resembling a new 
reconnection event (at t = 200sec) occurring in the 
scenario modified by the later ones. Table Q] shows 
the models that we use to discuss our results. We 
employed open and closed conditions for the lat- 
eral boundaries, whereas, initially, only a fix re- 
bound condition was used at the upper radial direc- 
tion, to simulate the outflows resembling the action 
of the reconnection site, and transmission condition 
for the bottom boundary, assuming that the pertur- 
bations are absorbed in the sunwards direction. Fi- 
nally, we also used open boundary conditions for 
the upper radial direction with a non-homogeneous 
initial background density (n = nosin(7ix/l) + m, 
(no,n{) = (2.07,2.53) 10 9 cm" 3 , x in Mm) to emu- 
lates the passage of previous SARO structures and 
to explore other plasma configurations. The charac- 
teristic boundary parameters were chosen in accor- 
da nce with typical obser ved dark lane structures, as 
in 
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Verwichte et al.l (2005) and iMcKenzie and Hudson! 
999). 



3. Results and Discussion 

3.1. SAROs plasma parameters 

Figure Q] shows, the Ml evolution of the three 
voided outflows at t = 196sec, it simulates, respec- 
tively, the voided density tracks left by three simulta- 
neous reconnection bursts, the temperature, the mag- 
netic field, the /^-parameter and the total pressure for 
open boundary con ditions. From Fig. [Tk we see that, 
as expected (e.g. Mc Kenzie and Hudsonl 1 199^) the 
vacuum density values are less than the environment 
ones in more than one order of magnitude (see the fea- 



Model 


B[G] 


n * 10 9 [cm 3 ] 


(AP/P)4 


Ml 


3 


0.46 


400 


Ml 


10 


0.46 


400 


M3 


10 


3 


400 


M4 


3 


0.46 


200 



Table 1: Simulated models: B the background mag- 
netic field in the radial direction, p the background 
density. 



tures above y = 20Mm); the temperature (Fig. [TJ)) is 
more than one order of magni tude higher in the vo ided 
cavities than outside (as in lLongcope et al. (2009) and 
I Warren et al.l (| 19991) ): the inner magnetic field is lower 
than its external value (< 1G) (Fig.[T]:), thus, the inner 
magnetic pressure cannot be responsible for prevent- 
ing the collapse of the vacuum zone. The /^-parameter 
(Fig. [TJl) is larger inside than outside. Furthermore, 
from the upper part of Fig. QJ we can see that the total 
pressure of the region is almost constant, i.e., the fea- 
tures that characterize the vacuum outflows have dis- 
appeared, meaning that the contour of the lanes are, on 
average, in total pressure equilibrium, as in Paper 3. 

The Fig. \T\ features are consistent with the dy nam- 
ical features described by Verwichte et al 1 (120051) ex- 
cept for the lack of oscillations in the 2D pattern and 
the homogeneity of the background medium. The 
multiple outflow dynamics reproduces the individual 
dynamic; i.e., in accordance with observational data, 
the temporal evolution analysis gave us decelerating 
speeds of the order of 100 - 300km s" 1 . The freezing- 
in of the plasma to the magnetic field induces the colli- 
mation of most of the kinetic energy towards the longi- 
tudinal direction, which is enhanced with the increase 
of the magnetic field. In accordance with Paper 3 (see 
Table 1 of the paper) -where we obtained a slow in- 
crease in the initial speed with the increase of the back- 
ground magnetic field- for a run of 200sec, varying the 
magnetic field from 3G to 10G (models Ml and Ml) 
the distance evacuated by the perturbation increases 
4Mm; while varying the background numerical den- 
sity from 0.46 1 9 cm" 3 to 3.0 10 9 cm" 3 (models Ml 
and M3) the distance evacuated by the perturbation di- 
minishes in 6Mm. 

At the initial stages of the flare event, before the 
formation of fan rays, the supra-arcade zone can be 
assumed as an almost homogeneous medium. FigureQ] 
intends to describe the initial dynamic of SAROs that 
descend from an apparent homogeneous background 
medium towards the arcades. Different background 
density give dynamical results consistent with the pa- 
ramete r range variations described by IVerwichte et al.l 
d2005l) . Thus, in the scenario proposed, the SARO dy- 
namics of several events is sustained by the interaction 
of localized (confined to the individual dark tracks) 
nonlinear waves and shocks acting in times compara- 
ble to the observations (see Paper 3). 

In the bottom region of Fig. [T] e.g., y ~ lOMm, 
we note an increase of the density (Fig. [T^), the J3- 
parameter (Fig. QJ1) and the total pressure (Fig. QJ) 
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Numerical density [10 9 cm" 3 ] 
0.02 0.36 0.70 1.05 1.39 0.56 



Temperature [MK] 
7.67 14.78 21.89 29.00 2.64 



Total Magnetic Field [G] 
2.73 2.83 2.93 3.02 0.0! 




Beta 
0.32 



0.44 0.55 0.35 



Total Pressure [dyn cm" 2 ] 
0.44 0.52 0.60 0.68 

■ I 



5 10 15 

Longitude [Mm] 



Fig. 1. — Simulation, at t = 196sec, of: a) density, b) temperature, c) total magnetic field d) ^-parameter and e) total 
pressure of the three outflows. 



with respect to their corresponding background val- 
ues; correspondingly we note a decrease of the mag- 
netic field (Fig. []]:). These bottom overdense features 
can be assoc iated with t he b right downward features 
described bv lMcKenzieU 2000). In our scenario, they 
are shock fronts that open the way towards the sun sur- 
face. They modify the background medium producing, 
among other effects, a mass pile up in front of the flow. 
A future 3D simulation could help to understand the 
relation and interaction of these reconnection outflows 
with the observational bright leading loops. 

The multiple shock dynamic structures produce 
temporal patterns w hich are reminiscen t of observa- 
tional constrictions dlnnes et al.1 12003 a|) . dlnnes et al.1 
l2003bh . see the knots in the SARO traces of Fig. [Q 
Figure |2] shows the total pressure at the heights indi- 
cated with cuts in Fig. QJ: ~ 32Mm, ~ 15Mm and 
~ 12Mm. Note that the total pressure is fairly constant 
in the upper region, near the triggering site at ~ 32Mm 
(dotted line). At ~ 15Mm, (solid line) the total pres- 
sure is higher inside the tracks. The constrictions or 
knots (e.g see the second outflow trace at ~ 12Mm, 
dashed-dotted line) correspond to regions of total in- 
ternal pressure lower than the outside values. 

As shown in Fig.[U in accordance with what is ex- 
pected from observations, it appears that the SARO 
lanes are not distorted by the fact that they are multiple. 




10 12 
Longitude [Mm] 

Fig. 2. — Total pressure of the three outflows at 
t = 196sec and heights: 32Mm (dotted line), 15Mm 
(solid line) and 12Mm (dashed-dotted line). See the 
constriction located at 12Mm in the second outflow 
trace. 
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There is not an appr eciable evidence o f the interaction 
between outflows. ISavage et al.l d2012l). in the frame 
of the scenario proposed by iLinton et all (|2009|) . in- 
terpreted the discreteness of outflows as an indication 
that the reconnection is highly lo calized. In our sce - 
nario, as in the one proposed by ILinton et al.l ( 2009 ). 
the production of plasma outflows is associated with 
reconnection bursts, but we do not assume a 2D config- 
uration nor a particular magnetic topo logy of the cur- 
rent sheet (Priest and Demoulin 1995). The scenario 
proposed bv ILinton et all (12009) relates inherently the 
plasma voids to reconnections that occur within the 
current sheet. In our case reconnection is not mod- 
eled but simulated as a pressure pulse that triggers the 
outflow production in an initially homogeneous media. 
Thus, we can still be concerned about the extent of the 
independent behavior between different outflows. 

3.2. Structure and wavy appearance 



An indication of the interaction between SAROs 
seems to be the wavy appearance that can be seen in 
some observations, e.g. IVerwichte et al.l (I2QQ5|) . There 
is an important difference between ID and 2D simula- 
tions in relation with the ability to generate the wavy 
appearance of the reconnection outflows. In Papers 1- 
2 (ID simulations) these features are simply obtained 
due to lateral rebounds that emulate a neighbor denser 
medium (see Fig. 2 and Fig. 1 of Paper 1 and 2 re- 
spectively). The oscillatory pattern is sustained over 
time and the wave structure does not lose energy. In 
Paper 3 (2D), the wavy appearance of a lateral cut is 
damped in few periods (see Fig. la of that paper), but 
if we focus in the 2D slides (see Fig. 2 of Paper 3) it is 
not possible to recognize the wavy character. It seems 
that the spatial distribution of energy of the 2D lateral 
rebounds (with closed boundary conditions) leads to a 
rapid damping of the oscillations and would not give 
account of the wavy appearance of the observations. 

Meanwhile, comparing the evolution of nonlinear 
waves produced by a triggering pressure pulse in the 
ID and 2D cases we recognize a similar dynamic be- 
havior that prevails during time and could explain the 
robustness of the voided patterns. Figured shows the 
spontaneous ID non-dimensional density evolution 
due to the action of a press ure pulse over a med ium 
with static initial condition dFernandez et al.lE009|) . In 
Figure [3k the triggering pulse is asymmetrically lo- 
cated, as in our runs. 

A detailed description of the evolution of the voided 



features can be obtained from these figures and from 
the wave analysis provided by the numerical tech- 
niques used. Figure [3k displays nondimensional time 
step numbers 200, 800, and 2000. Time step cor- 
responds to p(0) = 1 . In the figure we see two shock 
wave fronts for each time step, e.g., for step 800, the 
shocks are located at (x\ s - 39; x 2s = lll)Mm. Also, 
for the same time step two contact discontinuities are 
found at (x\ c = 46.5; X2 C = 103.5)Mm and two ex- 
pansion waves are found starting at (x\ e = 57; x 2e - 
93)Mm. Initially two shock fronts moving away from 
the triggering pressure pulse location are formed (see 
e.g. time step 200 in Fig. [3k and Fig. [3k). As the 
shocks travel along the magnetic field direction, the 
temperature and the density are increased making the 
energy rise, e.g. (x\ s = 39; x 2s - lll)Mm in time step 
800 of Fig. [3k- The density is abruptly diminished (e.g. 
(x\ c = 46.5; X2 C - 103.5)Mmin time step 800) by con- 
tact discontinuities that go behind the shocks and can 
be recognized because the pressure and the velocity of 
the flow are not changed while the waves pass. The 
temperature increases to maintain a constant pressure 
across the contact discontinuities. Also, two expan- 
sion waves -that are recognized because they dimin- 
ish the density while they pass- are initially produced 
and travel towards the triggering position, contrary to 
the shock fronts (e.g. {x\ e - 55.5; x 2e - 94.5)Mm 
at time step 200 in Fig. [3k). These waves collide, re- 
bound at the triggering location (see Fig. [3k-b, cen- 
ter part for time steps 800), and then travel in the 
opposite direction towards the contact discontinuities 
(e.g. (x ie = 57; x 2e = 93)Mm of time step 800 
in Fig. [3^). From the comparison of time step 800 
and time step 2000, we see that the expansion waves 
have a larger speed than the contact discontinuity. At 
time step 2000 the position of the shock fronts in 
Fig. [3^ are (x\ s = 7.5; x 2s = 142.5)Mm, the expansion 
waves have interacted with the contact discontinuities 
and lowered the density and pressure leaving behind 
a coupled nonlinear system of waves, e.g., the ranges 
(Ajci =31.5- 7.5; Ax 2 = 142.5 - 118.5)Mm in Fig. [3k 
and (Ajci = 28.5 - 15.0; Ax 2 = 109.5 - 91.5)Mm in 
Fig. [3k. As a result of this dynamic a central voided 
cavity is formed. 

Figure [3k, shows a 2D density cut taken along the 
SARO located at x = 19Mm in Fig.Hk. The similarity 
between Fig. [3k-b and Fig. [3k (compare t = 8 sec in 
Fig. [3k with time step 800 in Fig.[3k-b and t = 180sec 
in Fig. [3k with time step 2000 in Fig.[3k-b) expresses 
the likeness between the dynamic behavior of the ID 
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Fig. 3. — Upper panel) ID density evolution after 
a central triggering pressure pulse for different time 
steps, Middle panel) ID density evolution after an 
asymmetric triggering pressure pulse for different time 
steps, Lower panel) 2D SARO density evolution for a 
cut at x = 19Mm in Fig.[T^. 



and 2D descriptions. There are, however, small differ- 
ences in the shape of Fig. [3fc at t - 180sec from the 
ones in Fig.[3^-b, time step 2000. Comparing the fea- 
tures ranging from x = (10-20)Mm in Fig. [3fc with the 
corresponding features in Fig.[3k-b, we note than even 
when the patterns are similar the more irregular shape 
in Fig. 12: can be attributed to perturbations produced 
by the other SAROs. 

However, we found cases where the interaction be- 
tween SAROs is evident and has the observational con- 
sequence thaHtre^roa^esJhe wavy character regis- 
tered by IVerwichte et al.l (I2005|) . Figure Hfr-b shows 
the SARO evolution for Ml after a new pressure pulse, 
triggered at t = 200sec, has developed a wavefront 
over the discrete tracks left by the three initial SAROs 
that have traveled decelerating a distance of ~ 16Mm 
(Fig. |4k). In the upper part of Fig. 01) at t = 300sec it 
appears that the initial SAROs have lost their discrete- 
ness. We have stated that when the magnetic field is 
uniform the energy produced by the pressure pulse is 
distributed in a non-homogeneous way due to the col- 
limation generated by the magnetic field. While the 
SAROs propagate further a new trace left by the new 
burst can be seen sunwards (Fig. Ht). Also, the new 
reconnection pulse injects energy most of which is dis- 
tributed in a non-uniform way between the wakes left 
by the previous SAROs i.e., the initially spherical pres- 
sure pulse (Fig.|4k) at 212sec is later collimated by the 
magnetic field configuration at 300sec (Fig.|4j)). As in 
the other cases, the sunward elongation of the fourth 
lane is due to the interaction of upper rebounds and 
downward absorptions along the magnetic field that 
collimates the flow. 

A remarkable effect which has observational impli- 
cations is the distribution of kinetic energy. As can be 
seen from Fig. [5^, part of the new SARO kinetic en- 
ergy is distributed along the wakes left by the previous 
SAROs. In accordance with the observations, when 
the sequence is se en as a movie the wak e s exhibit a 
wavy motion (e.g. ( Verwichte et al.ll2005 ). dMcKenzie 



200(|, dKhan et al 



2007)). This can be explained due 



to the tendency of the wavefront energy to deviate into 
lower density regions, i.e., the wakes left by the pre- 
vious SAROs. The continuity Rankine-Hugoniot re- 
lation for hydrodynamic shocks, i.e., the shock in the 
longitudinal direction (Paper 1-3), and in the shock 
wave frame is 



Pd 

V u = — Vd 
Pu 



u upstream, d downstream. 



Initially, the energy of the pressure pulse is almost 
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spherically distributed due to its larger pressure value 
with respect to the ambient one (Fig. Hk), so we can 
consider that pd and Vd are almost the same all around 
the spherical boundary of the initial shock front. How- 
ever, the corresponding upstream p u values will be 
different depending on the previous SARO distribu- 
tion, i.e., p u is almost an order of magnitude lower 
in the wakes left by the previous SAROs. Thus, ini- 
tially the upstream flow speed will be larger in the 
SARO wakes than in the neighbor ambient. More- 
over, the previous SAROs have larger values of the 
flow speed than the ambient one and its value will 
be added to the upstream speed. This implies that 
two upstream shock front positions with initial values 
Pa < Pb, and thus va > v#, will continue their mo- 
tion with va > vb while the shock wave travels. Obvi- 
ously, the same relation stands for the kinetic energy, 
i.e. K A = Pav\ = Pbv b va > Pbv 2 b = K B (K the kinetic 
energy). Upstream, this relation is preserved as can be 
appreciated in Fig. [5^. Compare the kinetic energy of 
A and B in this figure with their corresponding densi- 
ties (pA < Pb) in Fig.[4t>. Note the correspondence be- 
tween the low-high density values in positions A and B 
in Fig. [5^ with respect to the high-low kinetic energy 
values in the same positions in Fig.0t>. The energy of 
the fourth reconnection event is partly deviated from 
the new wake into the previous ones depending on the 
background plasma parameter. Figure 0) shows the ki- 
netic density energy for M2 at t = 300sec. When the 
magnetic field intensity is increased (or the intensity of 
the pressure pulse is diminished, model M4) the wavy 
character is lost. 

3.3. A possible SAD configuration 

We performed a new simulation with a non- 
homogeneous background that emulates the modi- 
fication of the medium due to the action of several 
evolved SAROs. The previously simulated thin dark 
tracks (of sizes < 2Mm), emulate the action of re- 
connection pulses in a homogeneous medium. The 
bright overdense Shockwave fronts seen at the bottom 
of Fig. Q] would lead to a mass pile up in front of the 
flow changing the background density. To analyze 
the effect of new reconnections triggered in a medium 
modified by the SAROs, and that resembles the fan 
ray features described by the observations, we use a 
sinusoidal initial background density conditions. Fig- 
ure [6] shows the density variations obtained from the 
evolution of four pressure pulses located randomly in 
the medium distorted by previous SAROs. As shown 
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Fig. 4. — Ml numerical density after a fourth pressure 
pulse, the event is seen at left) t = 212sec, and at right) 
t = 300sec. Note the wavy features (zoom) due to the 
interaction between SAROs. 
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Fig. 5. — Kinetic density energy collimated into the 
wakes left by previous SAROs; left) Ml at t = 300sec, 
and right) M2 at t = 300sec. 
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in Fig. [5] the energy of new reconnection events is 
confined preferentially into the lanes left by previous 
reconnection events, i.e., thus the new dark voids are 
seen trailing brighter lanes (SAROs) once the medium 
is non-homogeneous. From Fig. [6] we note that the 
later pressure pulses, of a teardrop shape, are asso- 
ciated with different sizes (~ [2 - 9] Mm) depending 
on the effect produced by the new wavefront in the 
medium. Note the density range variation and the oc- 
currence of bright and dark lanes t hat resemble the fea- 
tures seen in Movie 2 provided by ISavage et al.l (120 12|) 
(see e.g. the 131 bandpass slice at 12 : 08 : 33UT). 
Also, the bottom slice at 12 : 22 : 09UT (Movie 2) 
exhibits a complex structure that resembles the action 
of evolved SAROs, where some loops are also rec- 
ognized at the left bottom. If this scenario describes 
a possible plasma configuration, SADs could be late 
reconnection events that occur in a non-homogeneous 
medium distorted by the action of the previous lo- 
calized SAROs and, the linking between SAROs and 
SADs would be related with the fact that energy of 
later reconnections is confined into the lanes left by 
the formers. This new burst of kinetic energy, suddenly 
confined into the older lanes, renews the accelerated 
motion of the SAROs and deforms their extremes until 
they acquire a teardrop-shape with an oscillating tail. 
Thus, these SADs would be seen trailing SAROs. This 
scenario does not exclude the possible presence and 
interaction between loops and SAROs as a result of 
the magnetic field reorganization, however, this needs 
of a 3D simulation that is beyond the scopes of this 
paper. A n important differenc e with the scenario pro- 
posed by ISavage etaD d2012|) is that all the obtained 
simulated voids, some of them resembling SADs, are 
not wakes but reconnection outflows formed by shocks 
and nonlinear interactions that could also be preceded 
by loops. In summary, our interpretation implies that 
SAROs and SADs are both reconnection outflows of 
the same type, excepting that SAROs are previous, 
and they modify the medium. Later reconnection out- 
flows will give the appearance that they are preceded 
by brighter lanes. What is lacking is to understand 
the interaction between SAROs and probable nearby 
loops. 

4. Conclusions 

We analyze the dynamic of multiple dark supra- 
arcade lanes that move decelerating towards the sun 
surface. In the scenario we proposed, dark lanes are 
explained as reconnection outflows (SAROs) result- 



Numerical density [10 9 cm" 3 ] 

0.09 1.37 2.65 3.93 5.20 

i i i 




5 10 15 

Longitude [Mm] 

Fig. 6. — New reconnection events occurring in a 
medium distorted by previous SAROs. 



ing from the evolution of shocks and expansion waves 
that form confined voided cavities -of high tempera- 
ture and ft values- collimated in the direction of the 
ambient magnetic field. These features sustain their 
morphology over time behaving as a whole pattern 
(almost independently from boundary conditions and 
other perturbations) with a sunward and decelerating 
motion triggered by upward reconnection events. 

When multiple SARO configurations are analyzed 
we note that the structure of each SARO is similar to 
the individual pattern reinforcing the non-interacting 
hypothesis. However, we found that the wavy charac- 
ter that can be seen in observations can be interpreted 
as an indication of interaction between SAROs. This 
interaction is significant when the bursts that trigger 
the phenomenon act over the lanes left by previous 
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SAROs. This wavy character is enhanced with the 
strength of the reconnection and/or with lower values 
of the magnetic field intensity. 

The fact that in a non-homogeneous medium i.e., 
modified by reconnection outflows (SAROs), the en- 
ergy of new reconnection events is collimated into the 
voided lanes left by the previous SAROs implies a 
linking between the previous and later reconnection 
events that can be consistent with the observations. 
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